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1 . SCOPE 


This doctiment describes the Yield Estimation Subsystem (YES) 
operational implementation of Robertson's Phonological Model 
for small grains. Specifically, this document describes the 
software established to compute Robertson's biometeorological 
time scale for all LACIE weather stations from approximations 
to daily maximum and minimum temperatures. This document 
does not describe the software at NMC that extracts the weather 
data or the software at JSC that writes the AI crop calendar 
reports . 



2 . APPLICABLE DOCUMENTATION 


1) Operating Procedures for the LACIE/YES Gridded Crop 
Calendar Report Writer by D.D. Wilcox. LEC-9955, 

Dec . f 1976. 

2) System Description of the LACIE/YES Gridded Crop 
Calendar Report Writer by D.D. Wilcox^ LEC»10004y 
Jan . , 1977 . 

3) AD 63-1347-4963^11 Documentation of Automated Crop 
Calendar Operation. 

4) "A Biometeoro logical Time Scale for a Cereal Crop 
Involving Day and Night Temperatures and Photoperiod". 
G.W. Robertson, Inti. J. Biomet ., 12#3 (191-223), (1968). 

5) Level 3 Baseline 

Yield Estimation Subsystem (YES) 

Requirements 

LACIE C00200 Sept., 1976 

6) Level 3 Baseline 

NOAA/NASA/USDA 

Interface Control Document 
LACIE C00710 Sept., 1976 

7) A FORTRAN Implementation of the Robertson Phenological 
Model, D.D. Wilcox, G. Champagne, R. Baskett , S . Wooley . 
LEC Tech. Memo, LEC-5974. 


A 



3. SYSTEM DESCRIPTION 


The Robertson Phenological Model is an algorithm for con- 
verting point estimates of daily maximum and minimum tempera- 
ture and day length to a point estimate of Robertson's Bio- 
meteorological Time Scale, a numerical scale describing the 
development of wheat. Knowledge of the stage of wheat develop- 
ment is useful to analyst interpreters (AI) in the estimation 
of wheat acreages. The Robertson Model has been implemented 
to make this knowledge available. 


The Robertson Scale ranges continuously from 1.0, called 
planting, to 6.0, called ripe. The Robertson Model assigns 
to each integral range 1.0 to 2.0,..., 5.0 to 6.0, an equa- 

tion fi,o-2.0^^max'^min' <^aylength) '^5. 0-6. 0 ^^max'^min 

day length) . To advance the model one day, let 

be the development yesterday. From ^ygsterday®®^®^^ 

Obtain yesterday's 


1.0-2. O' • 
^max'^min' 


*' ^5. 0-6. O' ^yesterday' 
and daylength, and compute: 


today ^^yesterday ^ ^yesterday ' 

Clearly t and t . are observations at the weather stations, 
max mxn 

It is most useful to have at the sample segments (LACIE 

modelling region, for example). Therefore, some method of ex- 
tending the weather station developments to the sample segments 
is necessary. 

It should now be unequivocal that any implementation of the 
Robertson Model must have three major components: 

1) A weather data collecting component, 

2) a crop modelling component, and 

3) a component that extends the output of (2) to the 
sample segments and writes a report. 

This document describes the structure and operation of com- 
ponent 2. Component 3 is described in documents 1 and 2 of 
section 2.0 of this report. Component 1 is in the implementa- 
tion stage. However, some software is described here which 
performs data collection and editing. It is anticipated these 
functions of component 2 will be abandoned when the final ver- 
sion of component 1 is put to use. From this point, component 2 
will be refered to as 'the system'. 



3 . 1 HARDWARE DESCRIPTION 


These programs and datasets are resident on the IBM 360/195 
complex at Suitland, Maryland. They should be transferable 
to any IBM 360-370 series machine with sufficient disk to 
handle the datasets and main memory to support the PL/I 
optomizing compiler. 


3.2 SOFTWARE DESCRIPTION 

The software consists of eight programs: one to prepare the 
station index file, six to prepare the crop calendar output, 
and one to copy this output from disk to tape. The program 
to prepare the station index file, DRECOVER, is input by 
cards at CCEA/Columbia. The six crop calendar preparation 
source programs are kept on cards at CCEA/Columbia and are 
listed in executable form on the partitioned dataset W.EDS. 
CCEA.PHASE3.LOAD; the executable programs are called by the 
following JCL files maintained on card decks at CCEA/Columbia: 

1) CLEAR, used to initialize disk file CROPDATA, 

2) CAPTURE, used to build the file, 

3) LIST, which lists CROPDATA and can be accessed at 
any time during preparation of the file. 

4) EDIT, which edits the file for mistakes and missing 
data, 

5) SUZYQ, which performs crop modelling, and 

6) INSERT, used to enter corrections into CROPDATA. 

JSCTAPE, the progrcun copying the crop calendar output to tape, 
is input by cards at CCEA/Columbia each time it is executed. 

In general, a crop calendar update will require CLEAR, followed 
by several runs of CAPTURE, followed by LIST. At this point, 
EDIT is run to correct for missing data and mistakes. SUZYQ 
is run producing updates; any errors not corrected by EDIT are 
detected by SUZYQ. INSERT is then used to complete editing 
of the weather dataset and SUZYQ is rerun. After the final 
run of SUZYQ, JSCTAPE is run to produce a tape for trans- 
mittal to JSC. 

Two (logical) files are maintained by SUZYQ to record the last 
predicted crop stage at each update, these files being in an 
Old Master/New Master relationship. In this way an error in 



the preparation of the data (ex. , incorrect data captured, 
not enough data captured) can be handled by going back to 
the Old Master and beginning the process anew. See Sec- 
tion 3. 2. 1.6 for further discussion. Data flow and file 
usage are shown in Figure 3-1. 
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Figure 3-1.— Data Flow and File Useage in Crop Calendar System 







Figure 3-1 (continued) 


DATASET NAMES AND FLOWCHART REFERENCES 


Flowchart Box # 
1 
2 

3 

4 

5 

6 

7 

8 
9 

10 


DSN 

CCEA.xxx.SRTD (see Section 3. 2. 1.1) 

CCEA. DAILY. XXX (see 3. 2. 1.2) 
W.EDS.CCEA.DAILY (see 3. 2. 1.3) 

W. EDS. CCEA. DAYDATA (see 3. 2. 1.4) 

W. EDS. CCEA. GROUP i.CROPDATA (see 3. 2. 1.7) 
W. EDS. CCEA. GROUPi. INDEX (see 3. 2. 1.5) 

W. EDS. CCEA. GROUPi.MASTERl (see 3. 2. 1.6) 
W. EDS. CCEA. GROUP i.MASTER2 (see 3. 2. 1.6) 
W. EDS. CCEA. GROUPi. GRIDINFO (see 3. 2. 1.8) 
W. EDS. CCEA. GROUP i. ARNO (see 3. 2. 1.9) 


XXX = month of interest, 3 letter abbreviation; ex. JAN, FEB, MAR 

i = I for Group I countries (U.S .A. , INDIA, CHINA) 

= II for Group II countries (U.S. S.R. , CANADA', 



The software components are divided into four classes: 


1) Weather Data Access Routines: 

These are a set of routines, INIT, SWAP, and TERM 
which provide a uniform method of reading and writing 
data to the weather data files. 

2) File Manipulation Routines; 

These routines include the main programs CLEAR, which 
initializes the weather data file, CAPTURE, which 
builds it, EDITOR, which edits out most mistakes and 
missing values, INSDEL, which makes card-input changes 
in the file, and PRINT which lists it. 

Also included is the program DRECOVER which sets up 
the INDEX file, the supporting routine CVB located in 
the PDS W. EDS. CCEA. PROD. LOAD and called by CAPTURE, 
and the JCL files which are needed for invocation of 
the main programs. 


3) 


4 ) 


The following is a list of the JCL files and their 
corresponding main programs ; 


File Name 

CLEAR ' 

CAPTURE 

EDIT 

INSERT 

LIST 


Program 

dLEAfe 

CAPTURE 

EDITOR 

INSDEL 

PRINT 


Crop Modelling Routine: 

This includes the main program CROPCALN, supporting 
routines START and PHENO, function DL, and the JCL 
file SUZYQ, necessary for invocation. 

Output Tape to JSC; 

This is a deck of JCL and the program JSCTAPE that 
produces the interface tape for shipment to JSC. 



3.2.1 TAPE AND DISK DATASETS AND ATTRIBUTES 

The following datasets will be referred to in subsequent 
descriptions of programs and subroutines. 

3 . 2 . 1 . 1 CCEA . XXX . SRTD 

North American data for past months on tape, where 'xxx' is 
the month abbreviation (ex., Jan. Feb., etc.) 

RECFM = FB UNIT = TAPE 9 

LRECL = 560 VOL = Eyyyyy, where 'yyyyy' 

BLKSIZE = 5600 is the tape number 

This dataset is accessed through the CAPTURE routine; a separ- 
ate tape exists at Suitland for each month of past data and is 
differentiated from the others by both the DSN and the tape 
number . 


3.2. ■'..2 CCEA. DAILY. xxx 

Non-North American data for past months on tape, where 'xxx' 
is the month abbreviation. 

LRECL = 384 
BLKSIZE = 12672 

All other information is identical to CCEA. xxx. SRTD discussed 
in 3. 2. 1.1. 


3. 2. 1.3 W. EDS. CCEA. DAILY 

North American temperature data for the current month on disk. 

RECFM = FB 
LRECL = 560 
BLKSIZE = 5600 

This dataset is accessed during CAPTURE. 


3. 2. 1.4 W. EDS. CCEA. DAYDATA 

Non-North American temperature data for the current month on 
disk . 


CRECL = 384 
BLKSIZE = 12672 



All Other information is identical to N.El^.CCBA. DAILY 
discussed in 3. 2.1. 3. 


3. 2. 1.5 W.EDS.CCEA.fGROJP I 'LiNDEX 

|GROO? 


RECFM > FB 
LRECL » 1024 
BLKSIZB « 10240 
DS0R6 » PS 

Yhe INDEX file for each group contains the number of stations 
on the file, followed by a list of station World Meteorological 
Organization (WHO) numbers and each station's three closest 
neighbor stations. 

INDEX is created by the program DRECOVER. 


3. 2. 1.6 W. EOS. CCEA.rGROUP I 1 ./master l) 

) GROUP II] |HASTER 2J 


RECFM * FB 
LRECL « 80 
BLKSIZE = 12960 
DSORG - PS 

Each crop calendar group must have two master files. The 
files are reversed each new two-week crop calendar period, 
such that during one two-week period MASTER 1 will be used 
as input and MASTER 2 as output; the next 14-day period will 
use MASTER 2 as input and MASTER 1 as output. This procedure 
is necessary to avoid loss of previously acquired data. 

The master used as input contains information on either the 
last day of the previous crop period or the winter wheat 
planting date if it has not yet been reached (winter wheat 
stations only) . 

If an error is encountered during execution of SUZYQ (ex., 
inverted or missing temperatures) , the mistake must be corrected 
and SUZJQ resubmitted, using the master files in the same 
oraer as the first run. 

MASTER 1 is initially created with an lEBGENER IBM utility 
program, using master cards on file at CCEA/Columbia as the 
SYSUTl input. 



3. 2. 1.7 


W.BDS.CCSA 


GROUP .CROPDATA 
GROUP II) 


RECFM ^ P 
LRECL = 1024 
BLKSIZE = 1024 
DSORG » PS 

This dataset contains the maximum and minimum daily air 
temperatures (or approximations thereto) for the 14 days 
for each station listed in INDEX. The CLEAR program is 
executed to change the starting day of the 14-day period^ 
and also to initialize all values to *9999'. The CAPTURE 
progr«un fills CROPDATA, and SUZYQ utilizes the file. 


3. 2. 1.8 W. EOS. CCEA. /group i\*GRIDINFO 

VgROUP II/ 

RECPM = FB 
LRECL = 80 
BLKSIZE = 12960 
DSORG = PS 

The GRIDINFO dataset is created or modified as output during 
execution of SUZYQ. The record contents are in the same for- 
mat as those of the output master file; while the output mas- 
ter file contains information on only the last day of a two- 
week period, GRIDINFO contains similar information for every 
day in the two-week period. GRIDINFO contains only information 
on stations with successful crop calendar executions; as a 
result, it is necessary to complete an error-free run of 
SUZYQ before this dataset is copied to tape using program 
JSCTAPE. 


3. 2. 1.9 W. EDS. CCEA. (“GROUP I > .ARNO 

|group iiJ 

RECFM = VBS 
LRECL = 796 
BLKSIZE = 800 
DSORG = PS 

This dataset is created by an lEBPTPCH step after SUZYQ is 
successfully run, and contains basically the Scune information 
as the output master file. It is accessed by JSC via a dial-up 
telephone line using the IBM Time Sharing Operation (TSO) . 



3.2.2 SUBROUTINE INIT 


INIT is called by CLEAR, CAPTURE, PRINT, EDITOR, INSDEL, 
and CROPCALN) INIT prepai.es the weather data file, CROPDATA, 
for direct access and initializes several variables used in 
SNAP. 


3. 2. 2.1 Linkages 

INIT defines COIMON blocks TALK, INDEX, DATA, and CLOSE, 
«diich are used by other subroutines and prograas. 


3 . 2 . 2 . 2 I nterfaces 

INIT must be called before SNAP. 


3. 2. 2. 3 Inputs 

INIT assumes that Unit 2 is attached to the INDEX file and 
that Unit 4 is attached to file CROPDATA. 


3. 2. 2. 4 Outputs 

After ejecting one page on the printer, INIT writes a heading 
for all output produced by subsequent subroutines and main 
program statements. 


3 . 2 . 2 . 5 Description 

INIT opens file CROPDATA for direct access, reads part of 
the INDEX file into the COMMON block INDEX and part into 
COMMON block CLOSE, prints a heading, creates CCX1MON block 
DATA, and initializes variables in COMMON block TALK, which 
is used in SWAP. 

3. 2. 2. 6 Flowchart 


Next page. 


3. 2. 2.7 Listing 
Follows flowchart. 











3.2.3 SHBROUTISE SWAP 


SWAP is the entry point used by all crop calendar prepara- 
tion routines to read from or vrrite to file CROPDATA. The 
'in' version of SWAP reads from CROPDATA. and the 'out' ver- 
sion writes to CROPDATA. 


3. 2. 3.1 Linkages 

SWAP uses the COMMON blocks TALK, INDEX, and DATA, which 
were declared by subroutine INIT. 


3. 2. 3. 2 Interfaces 


INIT must be called before SWAP. 


3. 2. 3. 3 Inputs and Outputs 


INIT assumes Unit 4 is attached to file CROPDATA. 


ACTION 


WMONUN 

TEMPS 

JULIAN 

RCODE 


if 1, then read from CROPDATA, with TEMPS and 
JULIAN as outputs to the main program; 
if 2, then write to CROPDATA, with TEMPS and 
JULIAN as inputs from the main program, 
station number of requested station, 
array of 14 days' maximum and minimum tempera- 
tures associated with WMONUM. 
first date of the 14-day period, 
if -1, then ACTION was other than 1 or 2; 
if 0, then transaction was successful; 
if 10, then WMONUM was not in the INDEX file. 


3. 2. 3. 4 Description 

SWAP maintains a BUFFER in COMMON block DATA to contain a 
16-station block of file CROPDATA; it also maintains a list 
of crop calendar stations in COMMON INDEX. The list of WMO 
numbers is in the same order as the data on CROPDATA. To 
access a WMO station, SWAP sequentially searches INDEX for 
the requested number. If the number is found, its position 
in INDEX is used to determine the block address (i.e., in 
which 16-station block of CROPDATA the station is located) . 

If the station number is not found, RCODE is set to 10. At 
this point it is determined whether the requested block is 
in core. If the block is in core, it is read from or written 



to; if the necessary block is not in core, it is swapped 
with the block currently in core. 


3. 2. 3. 5 Flowchart 
Next page. 


3. 2. 3. 6 Listing 


Follows flowchart 
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ORIGINAL PAGE 18 
OF POOR QUALITY 














RUN NO, 73 t 

JDE 5.C R I PT I ON 

MaeiTEii pile 
ftDQEO TO MftSTER 
LiST DATE COPIED 
LAST UPDATE 

PASSWORD 

PPO'jRAMMFR 

Language " 

PPOC parameter 


DATE 05/0A/77 TIME 1620 

SWAPS _ M ET DATA F OR AU TOMATION 

W. EOS. CCE A, GVM. master 
07/25/76 
> NONE 

07/25/76 — TTT9 

PDXS 

_ . VONHOLT 

For 

TNOJCL 


LISTING OF MODULE SWAP 


c**« 

1000 


1030 

1001 


SUBROUT I NE SWAP (ACTION , WMONUM « TFMPS 'ULl ANi I EK ) 
REAl*b BUFFER( 128 > 

INTEGER ACTION, WMONUM, TEMPSOA, 2 ) ,WM 0 ( 51 l) .POSIT 
^integer ACCFSS, CHANGE, page l.PAGFO 
INTEGER hlOCK.OFFsET 

iNTcGER"? NUvsTA,riLLEP.CCOATA( 2 . 16,161 
COMMON / T ALK/ incore, ACCESS -CHANGE, page I, PA 6 E 0 

common /inoe*/numsta,fillcr,wmo 

"COMMON” /OATA/CCDATA " " 

EQUIVAL£NCE(CCDATAI 1 ,1 , 1 ) ,flUFF£R(ln 

IFR s 0 

"IF' ( (ACTinMTUTnTT 0 R 7 (/TCTT 0 ls| 7 GT 72 n " 6 O 'TO "3000 

GO TO ( 1000 , 2000 ), ACTION 

SWAP IN 

00 800 J s 1 , NUMSTA 

IF UMO(J) "irOV WMONUM) — GO TO “900 

CONTINUE 
GO TO AOOO 
POSIT s J 

■ Block =■ (POSiT7r)“/"TS~~ r 

0 FFSET=P 0 SIT-(RL 0 CK- 1)*16 
IF (INCORE ,EO, RLOCK) GO TO 1050 
IF (CHANGE) 1040 , 1040 , 1030 

WRITE ( 4 « Income, 1001 ) BUFFER 

FORMAT ( 12 HA 6 ) 

PAGEO=PAGEO*l 

READ ( 4 *RLOCK, 1001 ) BUFFER 

'iNCOREsPLOCK 

PAGEI=PAGEI *1 

00 1060 J * 1 « 14 

TEMPS(JV 1 )sCC 0 ATA 117 J* 270 FFSETJ 

TEMPS ( J, 2) sCCOATA ( 2 ,J*2 .OFFSET ) 

JULlANsCCOATA(I,l, OFFSET) ; 

ACCESSsACCESS^l 





c«*» 

2000 


ACCESSsACCESS^l 

■PFTURN 

SWAP OUT 
DO 1600 J a 1 , 

IF (wMO(j) ,eo, 

‘CONTINUE 

GO TO 4000 
POSITsJ 



NUMSTA 

WMONUM) 


GO TO 


00000010 
00000020 
00000030 
_.00000040. 
00000050 
00000060 
00000070 
00000080 
00000090 
00000 00 
00000 10 
00000120 
"00000 30 
00000 .40 
00000150 
00000160 
“ 00000170' 
OOOOOIBO 
00000190 
_ P0000200 
00000210 
00000220 
00000230 
00000240 
00000250' 
00000260 
00000270 
_ 00000200 
00000290 
00000300 
00000310 
_ 00000320 
00000330' 
00000.340 
00000350 
_00000360 
00000370 
00000360 
00000390 
^00000400 
00000410 
00000420 


0T/2S/76 


C7/25/76 


RUN NO. 


DATE 05/04/77 


PLOCK = (POSIT-1) / 16 ♦ 1 

OFF5tT=POSIT-(RLOCK-l)»16 

IF (iNCOWe ,EO. MLOCK) GO TO 2050 
IF (CHANGE) 2040. 2040 * 2030 
2030 »WIT6(4 *inCORE.1001)8UFFER 

pageo=pageo*i 

2040 REAO(4«HLOCK.1001)9UFFER 
INC0RE=8L0CK 

PA G£ I = P A 6E_ 1 M 

'2050 OO 20o0 J = 1 . 14 

CCOATAd ,J*J?,0f FSET)=TEMPS( J.l) 

206 0 CCDATA(2. J»2. OFFSET )=TEMP S(J. 2) 
CCDATA (1.1. OFFSET) sJULI AM 
ACC£SS=ACC£SS*1 
CHANG£=1 

PFTjiPN 

'3000 *PI TE (6.3010)ACTION 
3010 format ( »0INVALID action C00E:».I7) 
CALL ERRTRA 

STOP 


4000 


lER = 1 
return 


C 


END 



3.2.4 SUBROUTINE TERM 


TERM is used by all crop calendar preparation routines to 
close the weather file and print statistics accumulated 
during the running of SWAP. 


3. 2. 4.1 Linkages 

TERM utilizes COMMON blocks TALK and DATA. 


3. 2. 4. 2 Interfaces 


SWAP must be called before TERM, i'^lRM is the last subrou- 
tine called before the end of any crop calendar preparation 
program. 


3. 2. 4. 3 Inputs 

TERM uses COMMON blocks TALK and DATA. 


3. 2.4. 4 Outputs 

TERM may write a block to CROPDAT A, assumed on Unit 4. The 
subroutine also prints out SWAP 'IN '/SWAP 'OUT' statistics 
and a closing message. 


3. 2. 4. 5 Description 

If a 16-station block remains in BUFFER, TERM writes the 
block to CROPDATA. Statistics accumulated during the execu- 
tion of SWAP and a closing message are written to the printer. 


3. 2. 4. 6 Flowchart 


Next page. 


3. 2.4.7 Listing 


Follows flowchart 



lucffe / 

XpA^viv^ S'tn.’fistity 
UflCLfcSi, 1 
\WAX.mx*. 





1.1 STING OF module TERM 


RUN NO, 73 DATE 05/0A/T7 TIME 1820 

DE«1C»IPTI0N TERM AUTOMAllON FILE. 


MasiTr.R TlLc. 

AOUEU TO MASTER 
LAST OiTc COPIED 
LAST UPDATE 

PAf.S*<ORO 

' programmer 
language 
p-<oc parameter 


W^EOS.CCE A, GVH. MASTER 
NONE _ 

none 

HLSK 

vonholt 

FOR 

tNOJCL 





SUBROUTINE TERM 

R£Aj,»t» 8 uFEER(128> 

INTEGER ACCESStCHANGE «PAGE 1 1 PA6E0 

COMMON /talk/ I NCORE , ACCESS . CH ANGE . PAGE I . PA6E0 

common /OATA/OuFFER 

IF (CHANGE) 200 « 200 t 100 
•*PlTt (A* INCORF.IIOBUFFFR 

format (]2BA8) " 

PAGEOsPAGEO*! 

WRITE (6.210 )ACCESS*PAGF. I *PAGEO 

format ( •-*• PAGING STATISTICS •*»./« ‘-ACCESSES* . 1 15./. 

'-PAGrTNS'TIlS./. *-PA6E OUTS' f'nA) — 

FORMAT('-*t/. »-•♦/«•-** END OF CAPTURE ROUTINE EXECUTION 
RE T URN 

"END 


' QOOOOQIQ 

„ 00000040 
OOOOOOSO 
00000060 
00000070 
OOOOOOSO 
OOOOOORO 

onoooloo 

OOOOQlIO 

00000120 

00000130- 

nnsui 





3.2.5 PROGRAM CLEAR 


CLEAR resets all entries in CROPDATA to aissing values (9999) 
and calculates the starting Julian date for the 14«day 
period about to be process^. 


3. 2. 5.1 Linkages 

CLEAR calls subroutines INIT» SWAP, AND TERM and uses COMfON 
block INDEX. 

3. 2. 5. 2 Interfaces 

CLEAR is the first crop calendar preparation routine to be 
run for any given 14-day period. 


3. 2. 5. 3 Inputs 

The INDEX file is on Unit 2 , CROPDATA is on Unit 4, and 
optional card input (date -reset value) on Unit 5. 


3.2. 5.4 Outputs 

The CROPDATA file is initialized. 


3. 2. 5. 5 Description 

CLEAR calls INIT to load COMMON block INDEX and open CROP- 
DATA to direct access. An att^pt is made to read a date- 
reset value from Unit 5; if no value is present, the reset 
value defaults to 14 days. The starting date of the 14-day 
period to be processed is reset and all temperatures in 
CROPDATA are set to *9999', using SWAP. The file is then 
closed by subroutine TERM. 


3. 2. 5, 6 Flowchart 


Next page. 


3. 2. 5. 7 Listing 
Follows flowchart 















3.2.6 JCL FILE CLEAR 


Executing JCL file CLEAR invokes the CLEAR prograa. 


3. 2. 6.1 Linkages 
None. 

3. 2. 6. 2 Interfaces 
None. 

3. 2. 6. 3 Inputs 

See Sections 3. 2. 5. 3 and 4. 2. 2.1. 


3. 2.6.4 Outputs 
None. 

3. 2. 6. 5 Description 

JCL file CLEAR binds the files referenced by prograun CLEAR 
to that progr2un and requests that the program be executed. 


3. 2.6.6 Listing 

//CLEARJtfPROC 

//CLEARbEXEC)iP(3«=CLEAR 

//STEPLIBJ4DDjiDSN=W . EDS . CCEA . PHASE3 . LOAD . DISP=SHR 
//FT02F001J<DDJiDSN=W.EDS .CCEA.f GROUPl"). INDEX, DISP=SHR 

1 GROOPig 

//FT04F001JiDDhDSN=W . EDS . CCEA . f GROOPI i. CROPDATA , DISP=SHR 

] GROUPIE 

//FT06F001JiDDJiSYSOUT=A * 

//FTO 5F0 0 lJ(SDD|tfDDNAKE=SYS IN 
//JtfhWWPEND 



3.2.7 PROGRAM CAPTURE 


CAPTURE reads the tape and disk tesqperature data files 
built at CCBA/Hashington and loads thesi into CROPDATA. 

3. 2. 7.1 Linkages 

CAPTURE calls subroutines IMIT, CVB, SHAP, and TEK^. 

3. 2. 7. 2 Interfaces 

CLEAR must be run before CAPTURE. 


3.2.7. 3 Inputs 

CAPTURE assumes the INDEX file is on Unit 2, the tape or 
disk data file on Unit 3, CROPDATA on Unit 4, and an optional 
card input file on Unit 5. This latter file is used with 
North American data only; it contains the Julian date of the 
first day of the month being processed, and the new year if 
the 2-week period spans tvo years. 


3. 2. 7. 4 Outputs 

CAPTURE loads CROPDATA with temperature data for the 14 -day 
period being processed. 


3.2.7. 5 Description 

CAPTURE reads the eleventh record from Unit 3 to determine 
if the data is "North American" or "non-North American". 

Unit 3 is rewound and INIT is called. For North American 
data, the Julian date of the first day of the month is read 
from Unit 5 and the array JULIAN is filled. North American 
data also requires that CAPTURE call subroutine CVB. For 
each new station SWAP (in) is called; CAPTURE determines what 
portion of the 14-day crop calendar is being accessed by this 
run of the program. CAPTURE makes the necessary adjustments 
if the 2-week period spans two years, including reading 
another card from Unit 5 for North American data to obtain 
the new year. Using SWAP (out), CAPTURE updates the CROPDATA 
file with the new weather data. After all stations have been 
processed, TERM is called and processing ceases. 



3.2.7. 6 Flowchart 


Next page. 

3. 2. 7. 7 Listing 
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3.2.8 JCL FILE CAPTURE 


Executing JCL file CAPTURE involves the CAPTURE program. 


3. 2.8.1 Linkages 
None. 

3.2.8.? Interfaces 
None. 

3. 2. 8. 3 Inputs 

Direct card inputs to this file are necessary only with North 
American data. See Sections 3. 2. 7, 3 and 4. 2. 2. 2. 


3. 2. 8. 4 Outputs 
None. 

3. 2. 8. 5 Description 

JCL file CAPTURE binds the files referenced by program CAPTURE 
to that program and requests invocation of program CAPTURE. 


3. 2. 8. 6 Listing 

//CAPTURE)z5PR0C)iUNIT= ,VOL= ,DSN= 'W.EDS . CCEA . DAYDATA ' 
//CAPTUREt4EXEC)6PGM=CAPTURE 

//STEPLIBj!5DDJiDSN=W . EDS . CCEA . PHASE3 . LOAD , DISP=SHR 

//FTO IFO 0 1)!5DD^SYS0UT= A , DCB= ( LRECL= 100, BLKS I ZE= 1 5 0 0 , RECFM=FB ) 

//FTO 2F00 l)!5DD)ziDSN=W . EDS . CCEA . f GROUPI 1L . INDEX , DISP=SHR 

UGROUPIlj 

//FT03F001Jj5DD}6DSN=&DSN. ,DISP=SHR,UNIT=&UNIT. ,VOL=SER=&VOL, 
//)ij!jj!}JzJLABEL= ( , , , IN) - . 

//FT04F001)i5DD^DSN=VJ. EDS. CCEA. VgROUPI \,CR0PDATA.DISP=SHR 

[GROUPI I T 

//FT06F001jzSDDjiSYSOUT=A 

//}6]6}6}6i6PEm 



3.2.9 SUBROUTINE CVB 

CVB is called by CAPTURE to convert WMO numbers from zoned 
decimal to binary (North American data onlj) . 


3. 2. 9.1 Linkages 
None. 


3. 2. 9. 2 Interfaces 
None. 


3.2.9._> Inputs and Outputs 

FROM the input WMO number (zoned decimal) 

NCHAR number of bytes in FROM (in this case, NCHAR=5 

TO the output WMO number (binary) 


3. 2.9.4 Description 

North American WMO numbers exist as zoned decimal and need 
to be converted to binary before processing can take place. 
CVB accomplishes this conversion. 


3 . 2 . 9 , 5 Flowchart 


Next page. 


3. 2. 9. 6 Listing 
Follows flowchart. 
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3.2.10 PROGRAM PRINT 


PRINT lists the contents of the CROPDATA file. It may be 
run at any tine. 


3.2.10.1 Linkages 

PRINT calls INIT, SNAP, and TERM and uses COmON block INDEX. 


3.2.10.2 Interfaces 
None. 

3.2.10.3 Inputs 

PRINT assunes the INDEX file is on Unit 2 and the CROPDATA 
file is on Unit 4. 


3.2.10.4 Outputs 

PRINT lists each of the stations on file INDEX and the 
corresponding temperature data from file CROPDATA. 


3.2.10.5 Description 

PRINT calls INIT and proceeds to call SWAP for each station 
on the INDEX file, printing the station number, starting 
Julian date of the current 2-week period, and the maximum 
and minimum temperatures in the 2-week period. After infor- 
mation for all stations on INDEX has been listed, TERM is 
called. 


3. 2. 1.6 Flowchart 
Next page. 


3.2.10.7 Listing 
Follows flowchart. 
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3.2.11 JCL FILE LIST 


Executing JCL file LIST invokes the PRINT progran. 


3.2.11.1 Linkages 
See Section 3.2.10.1. 


3.2.11.2 Interfaces 
None. 

3.2.11.3 Inputs 
None. 

3.2.11.4 Outputs 
None. 


3.2.11.5 Description 

JCL file LIST binds the files referenced by program PRINT 
to that program and requests execution of program PRINT. 


3.2.11.6 Listing 


//LIST0PROC 

//LISTJ!5EXECJ4PGM=PRINT 

//STEPLIB)[$DD)iDSN=lf . EOS .CCEA . PHASE3 . LOAD r DISP=SHR 
//FT02F001JtfDDj!iDSN=W.EDS.CCEA.CGROUPI I . INDEX,DISP=SHR 

(groupiiJ 

//FT04F001JiDDJtfDSN=W. EDS .CCEA. rcROUPI *1 .CROPDATA,DISP=SHR 

iGROUPIlJ 

//FTO 6F 0 0 1 JjDDliS YSOOT=A 
//tiWWPEflD 



3.2.12 PROGRAM EDITOR 


EDITOR checks for inverted t^peratures of missing data 
6uid attempts to replace bad values by substituting data 
from a given station's three closest neighbor stations. 


3.2.12.1 Linkages 

EDITOR calls, INIT, SNAP, CHECK, REPLAC, and TERN and uses 
COMMON blocks INDEX, VALl, and CLOSE. 


3.2.12.2 Interfaces 

EDITOR is not executed until all runs of CAPTURE have been 
ccxopleted and is run before the CROPCALN program is executed. 


3.2.12.3 Inputs 


EDITOR assumes that the INDEX file is on Unit 2, CROPDATA 
on Unit 4, and card input on Unit 5, containing the variables 
IP, IC, IE, and IS. 

These four variables attain the following values: 


IP = 0 no listing 

= 1 list the file CROPDATA 
IC = 0 do not call CHECK or REPLAC 

= 1 call CHECK but do not print results 
= 2 call CHECK and print results 
IE = 0 do not call REPLAC 
= 1 call REPLAC 

IS = 0 do not rewrite CROPDATA 
= 1 rewrite CROPDATA 


3.2.12.4 Outputs 


EDITOR produces corrections for CROPDATA and writes messages 
to the printer. 


3.2.12.5 Description 

Program EDITOR calls INIT and, using SWAP, fills arrays lAP 
and XTEMP. lAP contains the portion of the INDEX file con- 
tained in COMMON block CLOSE (the three closest neighbor 
stations for each station on the INDEX file) ; XTEMP contains 



file CROPDATA. EDITOR then reads a card, obtaining values 
for variables IP, IC, IE, and IS. Each tine a card is read, 
EDITOR incr^aents the variable IL and checks for end-of*file. 
The edit process is executed three tines. The first time, 
the contents of CROPDATA are listed; CHECK is called to iden- 
tify values that need to be replaced; REPLAC is called and 
messages are written to the printer. The second execution of 
the edit process calls both CHECK and REPLAC. During the 
third pass through the program, CHECK is called and CROPDATA 
is rewritten from XTEMP. TERM is called after the completion 
of the third execution. 


3.2,12.6 Flowhart 


Next page. 


3.2.12.7 Listing 


Follows flowchart. 
















Cck:'T 







•' • 0 | I w ; 7 > # I C-’ I V 

>l'l rJ~fU\j^^--^ FO-v Cko;-’ CALr'io.'-j EonisG 

-?FA|, XTt’*P<->il .JA.I)’. xmau<S 11.3) .ifSr>(Sil,31 
■~I']TrGrP'Tn'.»-T^irrj> .Hi;;-^TroTT73T%~riCTT51 i ) . 
I'lTF'oi-- n f^r ( i) ' •,*JIF '/ 

J J/ 1 / . )'J7/^/ ,.ni. I ^^(, J.J, IC. It . IS* (L 
I 7 ('^I 1 *.i) . |/wv(Sl I »i) 


I "j TFC't i ll )■ 1 s r A » F T LL ' 
l.7f>lCi\l «l TM.' S(SU * 1S.?> 

CO'iMiiij /IM'JEx/ 'lUMSTA.r ILLE^*/mO 
_C O M MOW / V a l_ I / X T K- 1 H ♦ X H A t# . X 50 . LO > 1 /•^ « H I r, Hp . 'Jj 1 1 AP* JUL T A M* TPLA C 

cn^MO-r /Y.r.ys^ /-^TAfM" ' ■ -■ * 

r\Li. iviT 

00 -lO t« = i*oir*STA 

CALL ( li-i*' lO'x.) *TE‘-'*3S. JiJLI AN) _ 

t.o 40 I = l.j 

TA*>(i<.i )=n 

no 10 j = 1 .NHMSTA 

_ I r ( V MO ( J) . h ) .31 ATM (K.n) GO TO <*** 

) CONTI -ait 
c»oro 3ii 
» f . I ) =.J 

) rp:. r J .yii ^ 

• /i *» 0 ' j’‘’= TTIT* 

00 t*V) I = \ 

IF (TFmPSIJ.I ) ,E 0,9v9) TEM^S ( j, I ) =99v9 

J } = T.f>?£S { J « I .) 

) ConTT^.’Ur. 

) r-r.’’i:-ut 

II =0 

- .*.< 7 . cr;'- a lu L _L£jil£L» H fJlS 

1 r O^ 1« I (ifii I <♦) 

!L=!L*1 

IF ( IL.OT.7) OTO 9000 


4^{fTATriVHr' 

FOi^mAII'I C"")- 

. • 'i .-H, >» r I . 'll r 

F^y^MA I < » « ) 

00 70 1 s 1.2 
j*i-l 

TTt (^ *1000.0) 

F Ok 1/4 I ( i « . 1 n « I 

ro ;TI .hi 
C O' 1 f J 'Or 


CAIAODI-9 OATA LISTEP — PASS*. 13) 
IM'^STA 


* A • I . 


Ji-Ji—I- t.lL.U.T^3) 

(\ 



(i )IO 

~ihi) 


rf.MPiK, 


P91 Tfc (0.1000'.) II 

F7.^'HATM CPn» CALANUK9 CMfeCPFR — PASS* .13) 

CAM OrCFMC.U) 

~^ F ‘-t i 7 - c ' r:i I '' vA t :S Oho 

• f 1 : (-, lf)f'.i->) |i 

M‘ C >'>••» CALAN0F9 FoiTOk — PASS*. 13) 
j: -'.LL ^-L/. C ( It j 

IF I I''.N‘' . I ) (MT'1 AO 

•.'ivl It (A. I000.>) Ti 

fowmaTi* Cpop Cai.AwOEP PE.VPITE FRF. — PASS*. 1 3) 

_0 0__0 2 1 1 _I< s I , ^ !I I V <^T /. 

C) ill) F = 
lO Jl'l J = 1 . ) A 

"T . 1 * . ) ^ 0 -. 3 ' "" ■ I ■ ■ 

cowTimjt 




CALL '5rtA>J(()iJT tTt JULIAN) 

C vJITt (*>»lfni-»7) ■.vMO(t') 

ir.0O7 FOWMAT(Mj»,I 6,» Rh ^WITTEN TO FILF.») 
3^‘l CONTINUE 
Noro ^^(> 

"•O'.tO CM I 7f 
ST'*^ 

t-. • If ■ 





3.2.13 JCL FILE EDIT 


Executing JCL file EDIT invokes the EDITOR program. 


3.2.13.1 Linkages 
None. 

3.2.13.2 Interfaces 


EDIT should not be executed until CAPTURE has successfully 
filled file CROPDATA. 


3.2.13.3 Inputs 

Direct card input to this file is necessary. See Section 

4. 2. 2. 3 for further discussion and input format. 


3.2.13.4 Outputs 
None. 

3.2.13.5 Description 

PROC EDIT binds the files referenced by program EDITOR to 
that program and requests the invocation of EDITOR. 


3.2.13.6 Listing 
//EDITJiPROC 

//EDIT)iEXEC)iPGM=EDITOR,TIME= ( , 15) 

//STEPLIB)iDDJiDSN=W . EDS . CCEA . PHASES . LOAD , DISP-SHR 
/ /FT02F001>5DD)iDSN=W . EDS . CCEA . fGROUPI . INDEX ,DISP=SHR 

1 GROUP I I J 

//FT04F001>SDD)iDSN=W. EDS .CCEA. fGROUPI 7 .CROPDATA, DISP=SHR 

(GROUPIlJ 

//FTO 5F0 0 1)!JDD)SDUMMY 
//FTO 6F 0 0 IJiJDDjiS YSOUT= A 
//}6}i)6]6}6PEm 



3.2.14 SUBROUTINE CHECK 


CHECK determines which data in CROPDATA should be replaced. 


3.2.14.1 Linkages 

CHECK uses COMMON block INDEX and VALl. 


3.2.14.2 Interfaces 
None. 


3.2.14.3 Inputs 


IC = 1 do not print results; 

= 2 print results 

IL the number of times through the edit procedure for 
this station. 

XTEMP (I , 14 , 3) an array containing the information on 

file CROPDATA, as well as space for the 
daily difference.' between maximum and 
minimtim temperature, where I is the 
number of crop calendar stations. 

IAP(I,3) an array containing the positions in 

array WMO of the three closest neighbor 
stations for each crop calendar station, 
where I is the number of stations. 


3.2.14.4 Outputs 


Outputs from CHECK contained in COMMON block VALl and where 
I is the number of stations invo^'/ed! 


XBAR (1,3) 
XSD(I,3) 
L0WB(I,3) 
HIGHB(I,3) 

NT (1,3) 
TFLAG(I,15,2) 


mean of each station's max, min temperatures 
and their differences. 

standard deviations of each station's three 
temperature types. 

lower confidence bounds of each station's 
three temperature types. 

upper confidence bounds of each station's 
three temperature types. 

number of acceptable values of each station's 
three temperature types. 

a logical array indicating "good" or "bad" 
temperatures; the 15th position i.s flagged 
"good" if all values for that temp, type are 
"good" and "bad" if at least one temp, is "had 



CDAY(14,2) array for printed output indicating values 
that should be nanually checked and edited 


3.2.14.5 Description 

CHECK begins by counting reason^U^le tai^ratares and flagging 
unreasonable ones. When a statxon has at least 5 "good" max 
t^ps. and/or 5 "good” min temps, for the t%#o>week period, the 
mean t^perature, variance, standard deviation, and upper and 
lower confidence bounds of the t^perature type are calculated. 
When both max and min a'' o present for a day, differences be- 
tween the two are calculated and summed; provided at least 
5 differences for a station exist, the same statistics as for 
the max and/or min are calculated. CHECK atteo^ts to find 
and calculate similar data for a station's three closest 
neighbor stations; 4-station averages of the mean temperature, 
standard deviation, lower and upper confidence bounds are 
figured. Whenever a large variation occurs between the mean 
of a temperature type and the corresponding 4-station mean, 
the 4 -station mean is substituted. 

Temperatures falling outside of the upper and lower confidence 
bounds are flagged. Then the mcucimum of one day is ccxopared 
to the minimum of the next; whenever too much or no variation 
between the two are encountered, temperatures are flagged 
according to the type of variation. After all t^peratures 
for a particular station have been checked, an array CDAY, 
denoting acceptable vs. not acceptable temperatures is filled. 
Beginning with the checks outside the upper and lower confi- 
dence bounds, the entire process is repeated. After the 
secord pass through this part of CHECK, the value of IC is 
checked. For IC=1, no printout is generated, but for IC=2, 
the following is printed: a listing of max and min tempera- 
tures and their differences, the number of values of each, 
the mean, standard deviation, and lower and upper confidence 
bounds of each temperature, as well as the array CDAY indica- 
ting values that should be manually checked and edited before 
procedure SUZYQ is run. 


3.2.14.6 Flowchart 


Next page. 


3.2.14.7 listing 


Follows flowchart. 
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3.2.15 SUBROOTINE REPLAC 


REPLAC att^pts to find and insert replac^ent values for 
t^pera tores flagged as unacceptable by CHECK. 


3.2.15.1 Linkages 

REPLAC uses COMMON blocks INDEX and VALl. 


3.2.15.2 Interfaces 
None. 

3.2.15.3 I nputs 

IE s 1 edit but do not print results 
= 2 edit and print results 

JULIAN ending date of the previous 2-week period 


where I is the number of stations involved. 


XTEMP(I,14,3) 

XBAR(I,3) 

LOWB(I,3) 

NT(I,3) 

See sections 3.2,14.3 and 
these variables. 


IAP(I,3) 

XSD(I,3) 

HIGHB(I,3) 

TFLAGd,] 0 , 2 ) 

3.2,14.4 for descriptions of 


3.2.15.4 Outputs 

Replacement values inserted into array XTEMP (1^14,3), where 
I is the number of stations, as well as messages indicating a 
successful substitution or the inability to make a substitution. 


3.2.15.5 Description 

This subroutine is executed twice each time it is called. 
REPLAC determines if any values for a particular station 
need to be replaced. When substitutions need to be made, 
several algorithms are available: 

A) At least seven temperatures of the station and type in 



question are available, as well as data from alter- 
nate stationer 

B) at least seven temps are available, but no alternate 
station dataf 

C) fewer than seven temps are available. 

Algorithm A will utilize up to four days of data from the 
primary station and up to five days of data from each alter- 
nate station, provided the stations also meet the criterion 
of at least seven "good" temperatures. For each station in- 
volved, the temperatures are standardized and averaged. 

A weighting factor is also calculated for each station, as 
well as an average weighting factor. The replacement value 
is the sum of the standardized value for the primary station 
and the average weighting factor. 

Algorithm B will first attempt to find a replacement value 
using the other temperature of the day if it is available 
(i.e., the max to find the min or vice versa). If this is 
not successful, REPLAC will try to use the average of up to 
four days of data for the replacement temperature. If no 
alternate days are available within a prescribed vicinity 
of the day in question, the replacement value will be the 
mean max or min value of the entire 14-day period. 

Algorithm C will use the average temperature of the alternate 
stations for the day in question as the replacement. If 
Algorithm C does not generate a replacement value, a message 
to that effect will be printed. 

After substitutions for all stations needing them have been 
made and the subroutine has been executed twice, control 
returns to program EDITOR. 


3.2.15.6 Flowchart 
Next page. 


3.2.15.7 Listing 


Follows flowchart. 


(.1 








iKJZ ff.GPtAC , CCA. T 


fVrT CaL| C*Vcwi«<\() ^6 

f#u CU^T^yit USiC 

*vlAc*-»vafc vrcxt»u^ tct^'>\ 
Stc^ <Vri%i t c« ,^rt »pcA ^4 1 »m) 
•♦yio* 4e»*«^5 *4*»rt*j 

0(c^ icic -fecttfr 


Cftcci;. {* 




n^rc 'i^A^icn) 







V? 

7 

re|5if^ceor>*^1 vetuiC 

*. rf^lA.« «..mrrf v«-W 

r -ffU-tcr 

JS 


_5 

Z 

nc^fac«^ uritu 

fic^uxrvivitiW v*»»Wy 
Ow.A0h(cf 
it 


g^/n- rvtv4 


SftA4 «t|T> 


XEx 2- 


lu*rr4<^ SubV't^fn 
\#i%^SS»uyL it 
















«5>H»0UT{«<P JlCUfl 

T^T^5|S^iM'^f^fir«Tr»imTrrniT**s^ 

iitrfH- jutiAN* iHL!.*|i: . _ ^ 

Of! *lii)0 l<»f« ^ * 

tf MjKiat,T#4 a-vi f HTA«is.i » » .*Nn, u»wT,TM.j>Qi,i$ aoTO #»«e 

f .l .M m w 1 nun r r . ... . 

ir i,^j0f,|ft*6(l«:TA,|0iiy«lM} I ««ttO 3900 
fr ♦NT4|ST*.tM|,tT,7> r.OTO 2200 


* CnftCi^ !'«•» S,t«T|ONS 



I » 1.3 

r** nA»rHrfl.ii,LT«ji t»o?o 20-<*o 

•'*-Tiri>tyr fy; T n ' « 'yrri 9^ 
l»f J‘4tt«.l»'i,l.T.7» OOTO 29*»d 
1 

; tH » t, < ■»{ tu 4 T t«) <HA»« 4 K 

20<»0 C(?"riu»4 *" * » • • 

V (H,U,U 0010 2U0 

^-wic r oi 'T' r ' ’SiiTW«^ ?r i 7 JrRigat g — 

70»/U/FI.O«T{.'0 


CHCCr fO» fctTC9*4*TC ST*T10N »I*3 
! J J4. 

7>l3iUu 
20ii>* j 

K' } if Sf f 


, ■■ - -n- - 

If c1rt4r^f*?ijiUYtiNH aoto 

* d ♦ ^ ^ f f f ,|MM 

MM 

tr Ctno nm% 

I ■*“ ‘ 1 ^* 11 ^. 

phhi iV**!NjLt. 2» onto 20»0 

2d 7.) ji|« y*.J0»V - — — 

2o*»0 jr>AV«j0A»*i 




I j^CfAv 14) .o». i 


?» jO . I • l.»f 
iraMAt till 

ir (TrL*4(i(»jo«y«iMH < 

d9^f> Cn‘itl»iufc 

im If ?n2,lt.^| §rt?5 M$5 


p4Y«!iiAt 




JIM MTO il» 


LT.4I iinfo tua 


)AV«|0At 

.mjtEj'.Mit 


2i%f iM,iT,ii lisa 



j2oa 7ftfe.a 


7? tPf!: jif lUtai •j'lTo 2210 
IP ( f f i*a I* M o*y • iW) i mw an 








» < I sTAJ • juLL iJ 


22 in 

c 

i "• 


'<»»»• ttfcfc'l.ACe MI$«>1NG V^tue WITH «5U»SITUTi VAI.UC 

- — ---.- 


70B?U**TfM*»(KfIOftY,lN4 

CONTlNUfc 

IF IN.L1.U SOTO 2€»00 
alternate station AVAItAHtE 
LOAT tN} 








3.2.16 PROGRAM CROPCALN 


CROPCALN reads the weather data and nodels the crop progress. 

3.2.16.1 Linkages 

CROPCALN calls subroutines INIT, SNAP, START, PHENO and TERM. 

3.2.16.2 Interfaces 

CROPCALN is run after file CROPDATA has been edited. 


3.2.16.3 Inputs 

The old master file is assumed on Unit 1, INDEX on Unit 2, 
and CROPDATA on Unit 4. 


3.2.16.4 Outputs 


Unit 6 is the printer, GRIDINFO is on Unit 8, the new master 
file is on Unit 9, and ARNO is on Unit 11. 


3.2.16.5 Description 

CROPCALN begins by writing a header to the printer and calling 
INIT. The program then enters a loop and reads a station. 

SWAP is called to obtain weather data from CROPDATA; CROPCALN 
checks and adjusts for a two-week period that spans two years 
and verifies the data obtained frcmi CROPDATA, writing error 
messages ceasing operations for that particular station when 
errors are encountered. If the crop stage has not yet reached 
1.0, START is called. If STAGE is between 1.0 and 6.0, 

CROPCALN goes into an inner loop; the daily crop progress for 
the two-week period is calculated via the subroutine PHENO 
and written to file GRIDINFO. Information is also written to 
file ARNO and the new roaster file. After all stations have 
been processed, ending messages are written to ARNO, the new 
roaster file, and the printer. TERM is called to print messages 
and close CROPDATA. 

3.2.16.6 Flowchart 


Next page. 


3.2.16.7 Listing 


Follows flowchart. 
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3.2.17 JCI FILE SUZYQ 

Executing JCL file SUZYQ invokes the crop calendar model 


3.2.17.1 Linkages 
None. 

3.2.17.2 Interfaces 
None. 

3.2.17.3 Inputs 
None. 


3.2.17.4 Outputs 
None. 

3.2.17.5 Description 

JCL file SUZYQ binds the files referenced by CROPCALN to that 
program and request the invocation of CROPCALN. 


3.2.17.6 Listing 

//SUZ YQ)!JPROCJiOLD= , NEW= , DSP=MOD 
//TRACKJ6EXECJiPGM=CR0PCALN 

//STEPLIB}^DD}tfDSN=W . EDS . CCEA . PHASE3 . LOAD , DISP=SHR 
//FT01F001|6DDJiDSN=W.EDS.CCEA.CGROUPI "I .MASTER&OLD . ,DISP=SHR, 

I^GROUPII f 

//J($J(5)iDCB=BUFNO= 1 

//FT02F001JtfDD}rfDSN=W. EDS. CCEA. /GROUPI *). INDEX, DISP=SHR,DCB=BUFN0=1 

USROUPIIJ 

//FT04F001)f5DDK5DSN=W.EDS.CCEA.5GROUPI .CROPDATA,DISP=SHR, 

I^GROUPIlj 

//j!5Jtfj!JDCB=BUFNO=l ^ 

//FT06F001bDDj!JSYSOUT=A 

//FT08F001bDDJtfDSN=W.EDS.CCEA./GROUPI .GRIDINFO,DISP=&DSP . , 

^GROUPIlJ 

//}^J6bDCB=BUFNO=l ^ 

//FT09F001)!JDDJiDSN=W.EDS.CCEA.i GROUPI 1 .MASTER&NEW. , DISP=SHR, 

[groupiiJ 

//}6WLCB=BUFm=l 

/ /FTl ir 0 0 lj!SDD)5DSN=W . EDS . CCEA . S GROUP 1 1 . ARNO , DISP=SHR , DCB=BUFNO= 1 

^GROUPIlJ 



//REPORT)rfEXECj!JPGM=IEBPTPCH . 

//SYSUTl)iDD)JDSN*W . EDS . CCEA . /GROUPI ( . ARNO , DISP=SHR 

IgroupiiJ 

//SYSIN)!$DD)5DSN«W.EDS .CCEA.DATALIB (PREFORM) ,DISP»SHR 

//SYSPRINTJiDD)rfDUMMY 

//SYSUT2)iDD)tfSYSOUT=A 

//)^)i)l$)^)[$PEND 
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3.2.18 SUBROUTINE START 

START acciunulates temperature data to determine the planting 
date for Individual spring wheat crop calendar stations. 


3.2.18.1 Linkages 
None. 

3.2.18.2 Interface s 
None. 


3.2.18,3 Inputs and Outputs 

CDEV cumulative development variable 
STAGE value of Robertson Scale 

CROP winter wheat = 11, spring wheat = 12 

MIN minimum temperature for a given day 

MAX maximum temperature for a given day 


3.2.18.4 Description 

Subroutine START automatically sets STAGE to 1.0 for winter 
wheat stations. Spring wheat stations must reach a CDEV>35.5 
before STAGE is set to 1.0. For average daily temperatures 
greater than 42°F, CDEV is increased by the value of 1.0. When 
average daily temperatures fall between 320F and 42°F, CDEV 
is increased by one tenth the difference between the day's 
average and 32°F. 


3.2.18.5 Flowchar t 
Next page. 


3,2,18.6 Listing 
Follows flowchart. 
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3.2.19 SUBROUTINE PHENO 


PHFNO is the routine which converts max and min te^>eratures 
and day length into crop progress via Robertson's model with 
Feyerherm's modifications. 


3.2.19.1 Linkages 

PHENO calls real function DL. 


3.2.19.2 Interfaces 
None . 


3.2.19.3 Inputs and Outputs 

TODAY Julian day being considered 

STAGE value of Robertson Scale 

CROP winter wheat =11, spring wheat = 12 

MIN minimum t^perature for a given day 

MAX maximum temperature for a given day 

MULT an array of 5 stage- and station-specific 

multipliers 

ALAT station latitude in hundreths of degrees 

ALONG station longitude in hundreths of degrees 


3.2.19.4 Description 

Subroutine PHENO calls function DL to compute the daylength 
of the day in question. The daylength. maximum and minimum 
temperatures are adjusted according to the stage of the pre- 
vious day. Three stage advance factors are calculated u.^ing 
stage-specific coefficients: VI uses daylength, V2 maximum 
cfei-perature, and V3 minimum temperature. AINCR, the day’s 
crop progress, is ccxnputed using the three stage advance 
factors and a station- and stage-specific multiplier; this 
is added to STAGE. If the crop stage is advancing too quickly 
during the fall. STAGE is held constant until a certain Julian 
date is reached. A check is made to see if the integer portion 
of STAGE has changed during the above calculations. If no 
change has occurred, operat'ons return to the main program. 

When change in the integer portion of STAGE does occur, the 
amount of chancn " ■'TT) is calculated, and the day’s progress 



is refigured. The value of STAGE to be returned to tliO main 
progreun is then calculated* adding the integer portion of 
STAGE to the product of the PCNT and AINCR. 


3.2.19.5 Flowchart 
Next i>age. 


3.2.19.6 Listing 


Follows flowchart 
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3.2.20 FUNCTION DL 

DL converts latitude and date to day length in hours. 


3.2.20.1 Linkage 
None. 

3.2.20.2 Interfaces 
None. 

3.2.20.3 Inputs 

XLAT station latitude in hundreths of degrees 
DATE Julian day being considered 


3.2.20.4 Outputs 
Day length in hours. 


3.2.20.5 Description 

After initializing constants, the function calculates the 
angle of the sun on the equator for the date and then adjusts 
the angle according to the station's latitude. The day length 
is computed using the product of the arc cosine of this angle 
and a constant value. 

3.2.20.6 Flowchart 


Next page. 


3.2.20.7 Listing 


Follows flowchart. 




^7 










3.2.21 PROGRAM INSDEL 


INSDEL is the main routine of the data replacement module. 
It allows the individual replacement of data in the weather 
file. 


3.2.21.1 Linkages 

Program INSDEL call subroutines INIT, SWAP and TERM. 


3.2.21.2 Interfaces 

INSDEL is run only the ^IT procedure cannot satisfactorily edit 
the entire contents of CROf^DAXA. If needed, INSDEL is executed before 
program JSCTAPE can be run. 


3.2.21.3 Ir^ts 

Ihe INDEX file is on Unit 2, CRDPDATA twi Unit 4 , and the card input 
containing data corrections is on Unit 5. 


3.2.21.4 Outputs 


Corrections are made in Q DEDATA. Sane error messages may ajpear. 


3.2.21.5 Description 

After calling INTT, INSDEL reads a card and scans it to see if it is an 
'INS'ert or a 'DEL'ete card. Using SWAP^ INSDEL then reads the span of 
days and type of tarperature to be changed. If 'DEL' v»as indicated, 
the appropriate tarperatures are set to missing (9999) ; if INS was 
called for, a card containing the new temperatures is read and thai 
tanperatures are inserted into CRDPDATA. Any time problans are encoun- 
tered with tlte card input, an error message is printed and the program 
proceeds to the next card. After all cards have been processed, TERM 
is called. 


3.2.21.6 Flovichart 


dext page. 

3.2.21.7 Listing 
Follows flowchart. 
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3.2.22 JCL m£ INSERT 


Executing JCL file INSERT invokes the IN93EL program. 

3.2.22.1 Linkages 
None. 

3.2.22.2 Interfaces 
None. 

3.2.22.3 Inputs 

Cards containing the infonnation necessary to correct JROIDATft oust be 
inserted. See Secticxi 4. 2. 2. 5 for formats. 


3.2.22.4 Outputs 


None. 


3.2.22.5 Description 

JCL file INSERT binds the files referaiced by INS3EL to that program 
£ind requests the executicxi of INSDEL. 


3.2.22.6 Listing 

//INSERTJdPRX 

//INSJ®ffiCJSPG»=INSDEL 

//STI3>LIBJd»)«KN=4tf .EDS .CCEA. HffiSE3 . IflAD ,DISP==SHR 
//FT02F001*JM»t6N=4(«.n>S.Ct3A.fGI^^ "L INDEX DISP=SHR 

(GROUPIlJ 

//FT04F001*DDW^DSN^.EDS.ClCEAyGR0UPI 1 .CKM1ATA,DISP===SHR 

IgrdupiiJ 

//FT05F0011tO)arOBtf1E===SYSIN 

//ET06F001*a»)zgdSYSOUT=A 



3.2.23 PROGRAM DRBOCVER 


This program fills the DOBC file from card iiQiut. 


3.2.23.1 Linlag es 
Ncx)e. 

3.2.23.2 Interfaces 


Whenever changes axe to be made in the HHEX file, DBEX3CWER is run to 
reconstruct the entire file. 


3.2.23.3 Inputs 

Input to DRBOOVER is by cards; the first card contains the nundaer of sta- 
tions to be in the file, punched in colxxois 1-3. subsequait 

card contains a staticxi to be in the file, and its three closest neigh- 
bor stat .OTs, in 4 (15, IX) format. 


3.2.23.4 Outputs 

The INDEX file is created. 


3.2.23.5 Description 

DRBCCfVER reads all input cards, sorting the data into NUM5TA (nutidaer of 
stations) and two cuxays: WO, containing all stations on the HDEX file, 
and STATN, containing the thr^ detest statiOTis to each station in 
array, VMO. DRBCXJVER then writes all data to the HCEX file. 

3.2.23.6 Flowchart 
Next page. 


3.2.23.7 Listing 


Follcws flowchart. 






//0»<EC0VER JOB ( '0010008E1HEA **»C 

//SF EXEC NCDOOCGO 
//SYS IN DO * 

Integer wMO(5ii),sTATN(5iit3) 


• » »COLUM« ) ♦V0NH0LT*REGI0Ns256K,TIME*1 


lNTEOtP*2 NUMSTA 
PEADISf l.EN0=2)NDMSTAt (WMO(I), (SrATI^(I*J) f J*l*3) 
FORMAT ( I3*/*2S5 (A < 15*1 X) */) *255 (A Hb* IX) ./) ) 

WRITE (8*3)NUMSTA* (WMOd ) * (STATNd* J) * J* I * 3) f l«l »5l IT" 
FORMAT(A2*2X*255AA*/,10(2d28AA) ./) ) 

STOP 

END 



3.2.24 PROGRAM J5CEAPE 


J5CTAPE is a PI/I program transferring data fron disk file GRIDIMX) to 
a tape to be s^t to Houston. 

3.2.24.1 Linkages 
None. 

3.2.24.2 Interfctces 

JSCTAPE shCHild not be executed until an error-free fun of has 

be^ processed. 

3.2.24.3 Inputs 

Disk file GRIDINEX) and a card containing header information for the 
t^)e. See Section 4.2.3 for card forma' 

3.2.24.4 Outputs 

JSCEAPE outputs a tape ccpy of file GRIDINFO preceded by a header indi- 
cating the countries and Julian dates involved. 


3.2.24.5 Description 

Tte program reads tlte card with header information and copies this infor- 
mation both to tape and to the printer. JSCEAPE proceeds to read each 
record of GRIDINFO and write the data to both the tape and the printer. 
When end-of-file is encountered, a closing message is written to the 
printer aid operation ceases. 


3.2.24.6 Flowchart 


Next page. 


3.2.24.7 Listing 
Follows flowchart. 
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4 . OPERATION 


Operation of the Robertson Phenological Model is a mixed 
huroan/machine procedure. This document describes the 
mechanical procedures. 


4.1 USER DOCUMENTATION 

Most of the work involved in operating the crop calendar is 
in the preparation of meteorological (met) data. To do this 
properly requires a knowledge of the organization of the crop 
calendar and met data files. 

Crop calendar countries are separated into two sets; Group I - 
the U.S., China, and India, and Group II - Canada and the U.S.S.R. 
Each set is run every two weeks with any given 14-day period 
ending on a Sunday, and the groups are run in alternating weeks 
(Group I one week. Group II the next) . Therefore the 14-day 
periods comprising each crop calendar run for each group are 
not identical. For example, the 2-week periods in Julian days 
for Group I may be 52-65, 66-79, and 80-93, while the 14-day 
periods for Group II may be 59-72, 73-86, and 87-100. 

Met data input is divided by months and by area (North America 
and non-North America). Previous months' data are on tape, 
while the current month's data are on disk. For each group 
there i- ; maximum of four met data collection runs per two- 
week periou. (last month, this month) x (North America, non- 
North America) - 


4.2 PROCEDURES 

In total, eight procedures are involved: one to construct 
the INDEX file, six to prepare data and run the Robertson 
Phenological Model, and one to prepare a tape of crop calendar 
data for use in Houston, 



4.2.1 CONSTRUCTING THE INDEX FILE (DRECOVER) 

Whenever changes are to be made in the INDEX file, DRECOVER 
is run to reconstruct the entire file. Card input follows the 
program cards as in the example below: 

//FT08F001)iDDJ6DSN=W.EDS .CCEA, fcROUPI \. INDEX, DISP=SHR 

IGROUPII) 

//SYSG0)6DD)6* ^ 

235 

42809^42971b42798Jd42492 

43003)i43192J643014Ji43117 

etc. 

/* 

// 


4.2.2 CROP CALENDAR PROCEDURES 

There are two sets of Job Control Language (JCL) PROC’s; 
one for the Group I countries and one for Group II. The JCL 
accesses the partitioned dataset W. EDS. CCEA.PHASE3 .LOAD con- 
taining all crop calendar programs in executable form. The 
deck for the appropriate group is input each time any crop 
calendar program is run. The EXEC card {or cards) for each 
program in the following sections of explanation refers to a 
PROC contained in these JCL decks, and is placed after the 
final //jiJJrfb^JJPEND card in the JCL deck. 

The JCL decks for Group I and Group II are found on the 
following two pages. 



JCL PROCEDURES FOR 


GROUP .1 



//LIST oPOC 

//LIST EXEC P6‘'=PPINT 

//STEPLIB DD DSf«!=K.EDS.CCEA.PHASE3.L0Ab*DISP=SHR 

//FTO ProOl DP 0«;M=te.Fr>c;.rrFft.-xoO!IPT:TNr)rx.nTS_P=c;HP 

//fJOUFDOl 00 0SN'=i(.E0S.CCEA.GR0UPI.CR0PDATA*01SP=SHR ^ 

y/FT05F001 OD SYSOUT=A 
// PEND 

//fott =>poc ^ ^ 

//EDIT EXEC P&m=E0IT0P*TImE=(,15) 

//STEPLT9 DD D5N=K. EDS. CCEA.PHASE3. LOAD. D1SP=SHR 
//FTOFFOOl DD DSN=W. EDS. CCEA. 6R0JPI . INDEX .DI SP=SHR 

//FT04F001 DD DSN=» .FPc;. CCF A. r.POtJPT . CROPPAT A . D7 SP=SMP 

//FT05FC01 DO DUMMY , • . - - 

//FT06FD01 OD SYSOUT=A - 

// DEND . - 

//SU7YQ PPQC 01 n=.ty;FW=.Dqp=Mnn ^ ^ ! 

//TPACK ExEC P6M=CR0PCALN 

//STEPLIB DD DSW=k'. EDS. CCEA.PHASE3. LOAD. DISP = SHR 

//FTOlFOOl DD DS*4 =W.EDS.CCEA.GPOUPI.mASTEPLOL 0. .DISP=SSP.DCe=aUFNO=l 
//FToprom r>n .rnc.rrFt..-,Dr>ttPi . «r,T «;p=cMP.r)rH=PiirK>r> =i 

//FT04F001 DD DSr>;=W.£DS.CC£A.5P0UPI.CR0P0ATA.DISP=SHP,DC9=BUFN0=l 
//FT06F001 DD SYS0J7=A 

//FT08F001 OD OSN=W.EOS. CCE A. 5P0UPI .6RIDI NFO.OI SP=LDSP. .DCB=BUFN0=1 

//FTQpFom nn nG\i=w.rnc;.rrF£.npnMPT -ma<;tfplnFj.^ .niGP=GHP.nrR=qtiFMQ=J 

//FTl JFOOl DD DSrv=lk.EOS.CCEA.GPOUPI.ARNO-DlSP=SHP.DC0=dUFNO=l 
//REPQPT EXEC PGM=IEBPTPCK 

//SYSUTI DD DSN=K. EDS. CCEA.GPOUPl. ARNO. D1SP=SHR 

//cystm on Ds\=k».Fns.rrf fl.naTfli TcupoFFoDM) .OT sp=SHP 

//SYSPRJNT DD DUMMY - , 

//SYSUT2 OD SYSOUT=A ~ 

// PEND 

_//CLEA^,T^ R PC:- ■■■: ^ ^ : 

//CLEAR EXEC PGM=CL£AR 

//STEPLIB DD DSN=W. EDS. CCEA.PHA5E3. LOAD. DISP=SHP 
//FTOEFOOl DD DSN=W. EOS. CCEA. GROUPI . INDEX »DISP=SHR 

on nsM=»'.'.t^ns.ccFA.GPniiPT .rooppATA .OT sp=SHP 

//FT06F001 DD 5YS0UT=A 
//FT05F001 DD DDNAME=SYSJN 
// PEND 

//I N'SEOT PRQC : : : ^ ^ I 

//If-'S E<EC PGM=INS0EL 

//STEPLIB DO D5N=W. EOS. CCEA. P-1ASE3. LOAD. D1 SP=SHR 
//FT02F001 DD DS'M= «<.•. EOS. CCE A . GROUP I . I NDE X . DI SP = SHR 

//F 7 Q4.rnm nn os\i=».i^ns.rrt ^ a .r.pouPT .roopnaTA . ni sp=smo 

//FT05F00] DD DDNAME=SYSIN 

//FT06F001 DD SYSOUT=A * ‘ 

// PEND 

_ /_/£ A P TURE PRQC t)NTT= . V01= . .i)SN= « W^FQ_S^Cf:FA.I)_A YOATA » ^ : 

//capture EXEC PGM=CAP7URE 

//STEPLIB DD DSN=W. EOS. CCEA. PMASE3. LOAD, DI5P=5HR 
//FTOlFOOl DD SYSOUT=A.DCB=(LRECL=100.aLK5I7E=1500.RECFM=FB) 

_ //f_I OPFnO I nn nSM=w.cnc; .rrrA . GPr>URT . iNprx .01 cptr^ rjp 

//FT03F001 DD D5N=LDSN. .OLSP=SHP.UNI 7=&UNl T. . VOL=SER=& VOL ♦ L ABEL= I . ♦ » 1 N ) ; 
//FTOAFOOl DD DS?v'=W. EOS. CCEA. GROUPI. CROPOATA.01SP=SHR 

//FT06F001 UD SYSOUT=A , 4 

yv P£ND ^ I i 










SVSOUT^A 


//FT05FOOJ 01 
// PEND 

//CLEAR PROC 
//n FAP FxEC 


//steplib do dsn=»«.eds-ccea. Phases. load, disp=5hr 

//riOEFOOl DD DSN=»k.E05-CCEA.GR0UPII.lNDEX.DISP = 5HR 
//FTOAFOOl DD DSN=W.EDS.CCEA.GR0UPI I .CR0PDflTA,DI5P=SBR 

//FTQSFOQI on 0nNAMF=SYSTN : 

//FT06F001 DD 5YS0UT=A 
// PEND 

//CAPTURE PROC UNIT=,VOL=-»DSN=«W.EOS.CCEA.DAYDATA* 

Z/C APJU?cL_i.AE C PGM =caPTUP E ^ : 

//STEPLlB DD DSN--W. EDS. CCEA.PHACE3. LOAD, DISP = 5HR 
//FTOlFOOl DD 5Y5OUT=A,0CB= (LRECL=100,BLK5I 7E = 1500,PECFy = FB) 

//F70PF001 DD D5N = w.E05.CCEA.&P0UPI I . I'MDEX ,D1 SP = SHR 

//FTnSr nni dp , pic;p = <;hp.unTT = ^ H NTT. .V01_ = SFR = X, VQ1 , LAP £L= 1 , . , J bJ 

//FTOaToOI DD DSIv=V;.EDS.CCEA.GR0U?II.CR0PDATA,01SP=SHR 
//FT06F001 DD SYSOUT=A 
// PEND 


4. 2. 2.1 INITIALIZING FILE CROPDATA (CLEAR) 


The CLEAR program is called by the CLEAR procedure. This 
sets all temperatures in W.EDS.CCEA.rCROUPI CROPDATA to 

\^GROUPIIJ 

'9999*. The starting date for the crop calendar period being 
run is also determined by CLEAR. The default option will add 
14 to the starting date of the last period to get the new 
starting date for the period being run. This may be over- 
ridden by us«r input, which is required if the previous two- 
week period spans two years. Examples are given below: 

A. The default option. 

//j^EXECjdCLEAR 

B. Setting the starting date 28 days after the last 
period. 

//)!5EXECbCLEAR 

//FT05F001)i5DDb* 

28 

/* 

C. Setting the starting date for a new year (in this 
example, the previous starting date is 362 and the 
next one is 010; the change is a -352 and the follow- 
ing input is used) . 

//)z5EXECbCLEAR 

//FT05F0Cl)i5DDK- 

-352 

/* 


4.2. 2.2 FILLING FILE CROPDATA (CAPTURE) 

The CAPTURE program is called by the CAPTURE procedure. This 
obtains maximum (MX) and minimum (MN) temperature data from 
catalogued disk files in the current month and from tape in 
past months. Four datasets can be accessed by CAPTURE: 

1. W. EDS .CCEA. DAILY - to access North American data in 
the current month (disk file) . 

2. CCEA. XXX . SRTD - to access North American data in a 
past month (tape access) where 'xxx' is the month 
abbreviation (ex., CCEA.JAN.SRTD) . 



3. H.EDS.CCEA.DAYDATA - for non-North American data in 
the current month (disk file) . 

4. CCEA. DAILY. XXX - for non-North American data in a 
past month (tape access), where 'xxx' is the month 
abbreviation (ex., CCEA. DAILY. JAN. ) . 

Data for past months is put on tape at Suitland. Maryland. 

The appropriate dataset name (and unit and tape numoer, if 
needed) is input on the EXEC card. When accessing North 
American data (W. EDS. CCEA. DAILY or CCEA. xxx. SRTD) the Julian 
date of the first day of the month that is being run must be 
input. This is not needed for non-North American data (W.EDS. 
CCEA.DAYDATA Or CCEA. DAILY. xxx ) . See examples below: 

A. GROUP I (China, India, and the U.S.) 

To CAPTURE data entirely within the current month, the 
following input is used: 

//JJEXECJrfCAPTURE , DSN= * W . EDS . CCEA . DAYDATA ' 

//)(JEXECJ6CAPTURE , DSN= ’ W . EDS , CCEA . DAILY * 

//FTOSFOOlbDDb* 

001 (eg., for January) 

/* 

To CAPTURE data whic’i spans a two*^nonth period, the follow- 
ing input is used: 

//bEXECJrfC APTURE , DSN= ' CCEA . JAN . SRTD ' , UNIT=TAPE9 , VOL=Exxxxx 
//FT05F0Cl)''ODbJt 
001 
/* 

//JrfEXECj«:APTURE,DSN= 'W.EDS .CCEA. DAILY ' 

//FT05F001)rfDDb* 

032 

/* 

//bEXECJrfCAPTURE , DSN = 'CCEA. DAILY .JAN ' ,UNIT=TAPE9 ,VOL=Exxxxx 
//bEXECbCAPTURE , DSN^ ' W . EDS . CCEA . DAYDATA ' 

/* 

2. GROUP II (U.S.S.R. and Canada) 

When running Group II the CAPTURE progreun that accesses 
W.EDS. CCEA. DAYDATA is on the CAPTURE procedure card in the 
JCL set-up ; therefore the EXEC card need only read: 

// J6EXECkJPGM=CAPTURE . See examples below. 

To CAPTURE data entirely within the current month, the 
following input is used: 



//)6EXEC)6>iPGM=CAPTURE 

//)6EXEC j!$CAPTURE , DSN= * W . EDS . CCEA . DAILY * 
//FTO5FOOl)llDD0* 

001 (ea., for January) 

/* 


To CAPTURE data which spans a twonnonth period, the 
following input is usedt 

//>Si:XECJSCAPTURE,DSN= 'CCEA. JAN. SRTD* ,UNIT=TAPE9 ,VOL=Exxxxx 
//FTOSFOOlJiDDJtf* 

001 (eg., for January as a past month) 

/* 

//>5EXEC)SCAPTURE , DSN= ' W . EDS . CCEA . DAILY • 

//FT05F0Cl)6DD)lS* 

032 (for February) 

/* 

//fcJEXECjiC APTURE , DSN= ' CCEA . DAILY . JAN ’ , UNIT=TAPE9 , V0L=EXXXXX 
//)iEXECjj5PGM=C APTURE 
/* 


C. Spanning two years 

Two runs of the CAPTURE program are needed; one will capture 
data for the end of the first year and the other will capture 
data for the beginning of the new year. 

When CAPTURE is used on North American datasets for the end of 
the first year, the Julian date tor December 1 is input as it 
was for previous December runs. Nothing else is required. 

When CAPTURE is used on North American datasets for the beginning 
of the new year, the Julian date for January 1, 001, is input 
as expected. On the next input card the value of the new 
year is punched; 

//)6EXEC)6CAPTURE ,DSN= ’ appropriate name ' 

//FT05F00I)!SDD)6* 

335 

/* 

//)6EXEC)iiCAPTURE , DSN= ’ appropriate name ' 

//FT05F001)!JDD)d* 

001 

197? 

/* 

When CAPTURE is used on non-North American datasets, no input 



is needed on either run because both the year and the Julian 
date can be read from the dataset. 

In all cases when two years are spanned, the year is needed 
to determine if the first year is a leap year so that the 
Julian dates oreater than 365 can be adjusted correctly. 


4. 2. 2. 3 EDITING FILE CROPDATA (EDIT) 

The EDITOR program is called by the procedure EDIT, and checks 
for inverted temperatures or missing data. There are three 
checks and two edits in this procedure. EDITOR scans the data 
twice and replaces or corrects missing data either by using 
data from the three closest stations or statistics generated 
oy the primary station to find a substitute. The third time 
EDITOR checks the data and prints it out indicating the sta- 
tions it could not correct, and stations where there are days 
with an unusually large or small gap between the MX and MN. 

The input for both GROUP I and GROUP II is as follows: 

//KJEXEOJEDIT 

//FTOSFOOlJtfDDJi* 

)6WO}SmiW)il)6WO 

Wi60l6WlW’MWl60 

WHOi 6 i 6 )i 2 i 6 i 6 ) 60 i 6 Wl 

/* 


The EDIT PROC needs 250K. 

To correct the data through input, use the INSERT procedure, 
explained in Section 4.2. 2. 5 


4. 2. 2. 4 RUNNING THE ROBERTSON PHENOLOGICAL MODEL (SUZYQ) 

SUZYQ is the Crop Calendar procedure. It executes program 
CROPCALN which uses the daily temperatures to compute how much 
the crop should have progressed each day for the 14-day period. 
SUZYQ prints out to three files: GRIDINFO, which contains the 
daily increments of the crop stage and is subsequently put onto 
tape; ARNO, which allows the summary data (the last day of the 
period) to be accessed from a remote terminal in Houston; and 
MASTER, which has two files, 1 and 2, that writes the daily 
increments to the "old" file and the stage of the 14th day to 
the "new" file. The input is as follows; 

//JiEXECbSUZ YQ , OLD= 1 , NEW=2 , DISP=OLD 
or //JiEXEC)tfSUZ YQ OLD= 2 , NEW=1 , DISP=OLD 



A record must be kept of which Master file (1 or 2) is used 
for input and output for each period (i.e.» for OLD and NEW). 

The I4th (or last) day of the period is written onto the "new" 
file eac** run, which on the next run becomes the "old" file. 

The Crop Calendar procedure will bcg^n incrementing at the 
correct date, using the 14th day of the "old" file to calcu- 
late the new starting date. The Master files are reversed 
each time SUZYQ is run. The SUZYQ procedure produces a sum- 
mary output which prints out the last day of the period ar ^ 
the crop's stage. 

It is easiest and most expedient to run EDIT and SUZYQ 
together. The EDIT "checker" Subroutine prints MX and MN, 
the difference between them, and indicates the days containing 
inverted temperatures or a wide or small gap between MX and 
MN by an 'XX*. SUZYQ prints the stations that have inverted 
temperatures and the day (s) on which this occurs (day 1, day 2, 
etc) in the 14-day period. These error messages appear in 
suimnary form as part of the crop calendar output. The input 
%40uld be: 

//bEXECkJEDIT 

//FT05F001bDD)i* 

wmnwiifiwii6wo 
wmwmnwiwtio 
Wl 60 i 6 W 2 ) 6 WOWtil 
/* ( 2 ) ( 1 ) 

//jiEXEOiSU Z YQ , OLD= 1 , NEW= 2 , D ISP=OLD 
/* 

When the two-week period spans two years, regardless of country, 
the value of the first year must be input: 

( 2 ) ( 1 ) 

//)tfEXEC)4SUZYQ,0i-D=l ,NEW=2 ,DISP=OLD 
//FTOSFOOlJrfDDJtf* 

1976 

/* 


As in the CAPTURE routine, the year is needed to determine 
if the first year is a leap year, so Julian dates in January 
can be adjusted correctly. 


4.2. 2. 5 MANUAL EDITING OF FILE CROPDATA (INSERT) 

If corrections need to be made for inverted temperatures or 
missing data, they are made with the INSERT procedure, which 
calls the program INSDEL. INSDEL will both insert and delete 
information; the two options are discussed separately. 


)I0 



A. To insert data 


//|l(EXEC]!fINSERT 

//PTOSFOOlJSDDji* 

INS|f$XXXXXj!l^t}000)((000 
aaaaabbbbli^cccddddd. . . . 

/* 

where: XXXXX is the station number 

mx,nm indicates whether the maximum or minimum 
temperature is to be changed 
000 the range of days, by Julian date, to be 

changed. If only one day is changed, the 
second field may be (xnitted. 
aaaaa.. the actual change, in fields of 5 (1415) 

An example — changing the maximum t^perature for station 
43758 to 61 and 63 on days 79 and 80 respectively: 

INSb4 3 7 5 SliMXbO 7 9k$0 a 0 
)6W6iai6i663 

If only one day for a station is changed, the input would be: 
INSb43758bMXJ5079 

Any number of INSERT change cards may be run at the saune 
time, in groups of two as shown in the above examples. 

When the INSERT procedure is run, SOZYQ should be run as 
a separate step after INSERT to be sure the changes were 
made and SUZYQ calculated the increments for those stations. 
One way to make a quick check is to make sure the last day 
of the period is the Scune Julian date printed out on the 
SUZYQ crop calendar summary output. The INSERT program 
alsc indicates if the correction was successful by messages 
routed to the printer. 

B. To delete data. 

//bEXECblNSERT 

//FTOSFOOlbDDb* 

DELbXXXXXJi^j^bO 0 ObO 00 
/* 

where: XXXXX is the station number 

mx,mn indicates whether maximum or minimum tempera- 
ture is to be deleted 

the range of days, by Julian date to be 
deleted. If only one day is deleted, the 

HI 


000 



second field may be omitted. 


An example: 

DELJtf 4 37 58)[iMN)[S08 3)[)08 5 
0ELJe(4 3758>$MXjfi083M085 

This example would replace the maximum and minimtun tempera- 
ture values on days 83 to 85 with *9999' (indicating miss- 
ing values) . 

As with the INS option, any number of delete cards may be 
run at the same time. 


4. 2. 2. 6 LISTING OF FILE CROPDATA (LIST) 

The LIST procedure calls program PRINT and can be executed 
at any time to list the contents of the CROPDATA file. In all 
the procedures that are run for the Crop Calendar (CLEAR, CAP- 
TURE, EDIT, SU2YQ, INSERT) the LIST procedure produces the 
listing of the data on CROPDATA and should be run with each 
program. 

The input would be: //Jt5EXEC)iLIST 


4.2.3 COPYING FILE CRIDINFO TO TAPE (DSJSCTAPE) 

This procedure, using program JSCTAPE, is run to put the crop 
calendar daily increments from file GRIDINFO onto tape to be 
sent to Houston. The JCL for DJSCTAPl (for GROUP I) or 
DJSCTAP2 for (GROUP II) follows; 

job card 

/7J6bJeJTIME=l , CLASS=C 
//bEXECJtJNPLIXCLG 
//PLI .SYSIN)6DDJlJ* 

//GO.INbDDbDSN=W.EDS.CCEA. CgROUPI ^ .GRIDINFO, DISP=SHR 

(groupii) 

//G0.0UT)iDDJrfUNIT=TAPE7,DISP=(NEW,KEEP) , LABEL= ( 1 , BLP , ,OUT) , 
//16}6}6DCB= ( RECFM=FB , LRECL= 8 0 , BLKS I Z E= 3 2 0 0 , DEN= 2 , TRTCH=ET ) , 
//)^bbVOL=SER=Exxxxx 
//GO.SYSINbDDb* 
label card 
/* 

On each run the VOL=SER= ' tape number' must be changed and the 
appropriate date for the period must be put on the label card 



following the card //GO.SYSINJ^DDjrf*. The format for this card 
is as follows: 


columns 1-6 
7-8 
9-11 
13-14 
15-17 
19-80 


CCCOUT 

year of first day of the run period 
Julian date of first day of run period 
year of last day of run period 
Julian date of last day of run period 
commentary - used to indicate the GROUP 
and countries on the tape 


example: CCCOUT77101)rf77114)iGROUPII-COUNTRIES :USSR, CANADA* 


After the DJSCTAP job has been successfully run, the tape 
must be sent to Houston. CCEA/Columbia calls Suitland and 
gives Mr. Pennington the tape number, the group number, and 
the Julian dates for that run period. This information is 
for the outside label. 


4.3 FLOWCHART 

A flowchart showing the sequence of the crop calendar proce- 
dures is found on the following page. 














